Dynamics of particle deposition on a disordered 
substrate: I. Near-Equilibrium behavior. 



£ 

O 

o 



X 



Yan-Chr Tsai [] and Yonathan Shapir 
Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627 



D 

00 : Abstract 



^1 . A growth model which describes the deposition of particles (or the growth 

of a rigid crystal) on a disordered substrate is investigated. The dynamic renor- 

ON 



malization group is applied to the stochastic growth equation using the Martin, 
Sigga, and Rose formalism. The periodic potential and the quenched disorder, 
upon averaging, are combined into a single term in the generating functional. 
Changing the temperature (or the inherent noise of the deposition process) two 



different regimes with a transition between them at T sr , are found: for T > T sr 
this term is irrelevant and the surface has the scaling properties of a surface 
growing on a flat substrate in the rough phase. The height-height correlations 
behave as C(L,r) ~ ln[L/(r/L 2 )]. While the linear response mobility is fi- 
nite in this phase it does vanish as (T — T sr ) 1,78 when T — > T+. For T < T sr 



^Present address: Department of Physics and Astronomy, University of Pennsylvania, Philadephia, 
PA 19104-6396 



there is a line of fixed-point for the coupling constant. The surface is super- 
rough: the equilibrium correlations function behave as (In L) 2 while their short 
time dependence is (lnr) 2 with a temperature dependent dynamic exponent 
z = 2[1 + 1.78(1 — T/T sr )\. 

While the linear response mobility vanishes on large length scales, its scale- 
dependence leads to a non-linear response. For a small applied force F the 
average velocity of the surface v behaves as v ~ F 1+< *. To first order ( = 
1.78(1 -T/T sr ). At the transition, v ~ F/(l + C\ ln(F)|) L78 and the crossover to 
the behavior to T < T sr is analyzed. These results also apply to two-dimensional 
vortex glasses with a parallel magnetic field. 
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I. INTRODUCTION 



Much progress has been achieved recently in the understanding of surface growth in processes 
of deposition, sedimentation, epitaxial growth, solidification, etc A few years earlier the 

static and the dynamics properties of roughening of crystalline surfaces were elucidated [f||j|. 
Recent investigations have concentrated on the connections between surface roughening due to 
thermal fluctuations on one hand and that due to the kinetic growth itself, on the other hand. 
The underlying discrete structure of the particles (or the lattice) may lead to a kinetic phase 
transition between smooth and rough phases or between two rough phases with distinct scaling 
properties [![]. Since the underlying discrete structure is relevant at low temperature (or a low 
noise regime in the deposition) one cannot escape the question on how disorder in the substrate 
might modify the surface properties. 

The effect of the substrate disorder on the dynamics of the growing surface is the subject of 
our analysis ||. We address this issue using the dynamic renormalization group (RG) applied 
to the stochastic growth equations using the Martin, Sigga, and Rose (MSR) formalism 0. 

In the present paper we address the dynamics near-equilibrium. This regime is character- 
ized by a very slow rate of deposition such that the system is very close to thermodynamic 
equilibrium. In this regime only slight modifications from equilibrium are considered. In par- 
ticular the fluctuation-dissipation theorem (FDT) [^,[| and the Einstein relation ( between the 
mobility, the diffusion constant and the temperature) both hold. 

Far-from equilibrium, the growth equation does not obey the FDT. The symmetry under 
h — > —h (h is the height of the surface) and time reversal t — > — t are broken. The most relevant 
additional term, as shown by Kardar, Parisi, and Zhang (KPZ) |10| is due to the lateral growth 
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of the oblique surface. The behavior far-from equilibrium will be the subject of a following 
paper. 

In general, the scaling properties of the growing surface are manifested in the height-height 
correlation function: 

C(L, r) =< (h(x + L,t + r)- h(x, t)f > (1) 
or the corresponding surface width: 

W(L,t) = [C(L,t)] 1 ", (2) 
which obeys an asymptotic behavior of the form: 

W(L,r) = L a f(r/L z ). (3) 

In this expression a is the roughness exponent which characterizes the extent of the rough- 
ness of the surface and z is the dynamic exponent. f(x) is a scaling function which approaches 
a constant for large x. For small x ( r << L z ): f(x) ~ x 13 where (3 = a/z. At early stage of 
the growth the surface roughness increases as W{t) ~ r 13 while for r >> L z , W depends only 
on L and behaves as L a . 

In the absence of any disorder in the substrate the near-equilibrium behavior was analyzed 
extensively in the context of surface roughening. The original work was due to Chui and Weeks 
(CW) 0] and this system was further analyzed by Nozieres and Gallet (NG) 0. Their most 
important findings were as follows: 

In the high temperature rough phase 

C(L,T)~\n[Lf(r/L z )] (4) 
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which corresponds to a = 0, (3 = and z = 2. The scaling form of the correlation function 
in Eq. @ can not apply to that of Eq. (f|), since both a and j3 vanish but come with a finite 
ratio. 

In this phase the effect of the discreteness (or the lattice) is not relevant. This behav- 
ior is equivalent to that of a free-surface in which the surface tension is the only interaction 
determining its properties. 

In this regime the macroscopic mobility defined as the ratio between the average velocity 
v =< ^ > and the "force" F driving the surface, is finite. 

In the smooth phase C(L) is independent of L, and the mobility vanishes. The mobility has 
a finite jump from a finite value to zero at the roughening temperature. The growth process 
at low temperatures is by nucleation of higher "islands" on top of the smooth surface. This 
"activated" growth has drastically different dynamic properties which are determined by the 
diffusion of the deposited particles on the surface and their attachment to the "islands". 

While we study here the surface properties in deposition of cubic (or tetragonal) rigid 
particles, our study applies as well to the growth of crystalline surfaces if the rigidity of the 
solid is large enough. Our theory will apply if the surface height is smaller than the scale on 
which the random deviations in the substrate cease to affect the positions, along the growth 
direction, of the lattice ions. This scale will be larger the stronger is the Young's modulus 
which measures the longitudinal rigidity. 

It turns out that the same stochastic equation of motion also describes the behavior of other 
random two dimensional systems. The most important case is a system of vortex lines in a 
superconducting film with the applied magnetic field parallel to the film. (Charge density waves 
at finite temperature is another such a system). Therefore the conclusions of our investigations 
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also apply to 2D vortex-glasses [p| -[T^|. We shall come back to these implications in the last 
section. 

The outline of this paper is as follows: In the section II we present the stochastic equation 
of growth and the related MSR generating functional. In section III the RG scheme is outlined 
and the recursion relations are derived. Section IV is devoted to the discussion of the results 
and their physical implications. The final section V is dedicated to a summary of the important 
conclusions. In the appendices we provide more details of the RG calculations. A short letter 
announcing the most important results was published elsewhere M. 

II. THE EQUATION OF MOTION AND THE ASSOCIATED GENERATING 

FUNCTIONAL 

The prototypical paradigm for the simplest deposition process is the Edwards and Wilkinson 



(EW) model [15| for the sedimentation of granular particles. The continuum limit form of their 



equation of motion for the height h(x, t) is: 

d ~ 
/T 1 — t) = vV 2 h{x, t) + C(x, t)+F (5) 

jl is the microscopic "mobility" of the upper surface, v is the " diffusion constant" for the 
particles on the surface, F is proportional to the averaged deposition rate which is very small 
(large deposition rate will be discussed in a second paper), and ((x,t) is the local fluctuation 
from the averaged deposition rate, which obeys: 



< ((x,t)((x,t) >= 2DS'\x-x)S(t -t). (6) 
We can define the effective " temperature" of this system by the Einstein relation: T = Djl. 



If the discrete nature of the particles is taken into account the height h of every column of 
particles must be an integer multiple of the vertical size of the particle b. This discrete constraint 
leads to a periodic 5-function potential on h. This periodic potential may be expanded in 
Fourier series of which only the basic harmonic is relevant. On a flat substrate there will be 
an additional term of the form 7^ sin (jh^, t)) on the r.h.s. of Eq. |5| (7 = ^ and is the 
amplitude of the periodic potential). 

In the presence of a random substrate the minima of the potential will be randomly, and 
independently, shifted for each column. Hence the equation of motion becomes: 

jl-^hix, t) = uV 2 h(x, t) - 7^ sin[7[/i(f, t) + d(x)\\ + ((x, t) + F (7) 

d(x) is the local deviation of the disordered substrate as depicted in Fig. 1. The associated 
"phase" Q(x) = 2nd ^ is uniformly distributed between and 2tc, and is uncorrelated for 
different location x on the substrate. 

Then the equation of the growth process becomes: 



ji -1 ^ = F + u(V 2 h(x,t)) +^sin[ 7 h(x,t) + 6(f)] + C(x,t), 



where a is the lattice constant in the horizontal plane. 

To investigate this stochastic equation systemically, one can utilize the MSR (Martin, Sigga, 
and Rose) formalism by introducing an auxiliary field h to force Eq. (H) through a functional 
integral representation of a 5-function. The generating functional for Eq. (|8|) takes the form as 
(after averaging over ((x, t)): 

Z e [J, J}= J VhVhexp{S [k h] + Si[h, h] + J d 2 xdt(Jh + Jh)}, (9) 

where 



S [h, h}= [ d 2 xdt[Djl 2 h 2 - h(—h - fiuV 2 h)), (10) 

Si = [ d 2 xdt[^-hsm{ 1 h{x,t) + 9(f))]. (11) 

The generating functional Z[J, J] can be directly averaged over the quenched disordered d{x) 
because Z[J = J = 0] = 1 [pOfl . One may calculate any averaged correlation and response 
function by differentiating the generating functional with respect to the current J or auxiliary 
current J and setting J = J = 0. 

After averaging over the disorder the effective generating functional reads: 

(Zq[J, J}) disorder = J VhVheX P { J d 2 xdt[Djl 2 h 2 - h( — h - jiuV 2 h)) + 
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d xdtdt'h(x, t)h(x, t') cos(7(/i(x, t) — h(x, t'))}, 



(12) 



where g = y 2 . If we choose D = T/2 -1 , the system will evolve into the configurations weighted 
by a Boltzmann factor e~ H ^ T , which obeys the fluctuation dissipation theorem (FDT) ||. To 
simplify the calculation, we redefine those physical parameters as: fxv = fi,Du~ 2 = D,g = 
gu 2 , v = uu 2 , F = — . Then the equation of motion becomes: 

x dh{x,t) =F+ 2 + 7| gin(7 ^ (f) t) + + c(f; £) _ (13) 

at a z 

Here (((x~i, ti)£(x~2, t^)) = 2D6^(x\ — x~2)6(ti —t 2 ). The resulting generating functional for the 
present case reads: 

Ze[J, J}= J V4>V0exp{S o [4>, 4>] + S^, 0] + / d 2 xdt(J4> + J (/>)}, (14) 

where, 



S [<f>,<t>] = / d 2 xdt\D^ 2 -^-/iVV)], (15) 



SV = / d 2 xdt[^0sin( 7 0(f,t) + 9(f)], (16) 

where h(x,t) = <fi(x,t) and h(x,t) ~ 4>(x,t). In the same way, we arrive at the averaged 
effective generating functional: 

~,d 



(Z e [J, J}) dis0 r d er = J VfrDQ exp{ J d 2 xdt[D /i 2 4> 2 -$(—<j>- /iV 2 0)] + 

d 2 xdtdt'(p( y x, t)(f){x, t') cos(7(0(f , t) — (p(x, t'))}. 

(17) 
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In the generating functional in Eq. (|TTD, the term which contains the cos[7(0(x, t) — 0(x, t')] 
is non-local in time and will be responsible for the creation of non-trivial Edwards-Anderson 
correlations. 

In the renormalization process it turns out that this term generates a new " quadratic" 
term, also non-local in time, of the form: 

If 

-/i 2 z/ J dxdtdtV<j)(x,t)V<j){x,t') (18) 

We therefore add this term to the generating functional and will follow the flow of v as well 
(as shown below it will play a crucial role in altering the long-range height-height correlations.) 

Here we treat the the last term in Eq. (|H]) as a perturbation to the free action, and expand 
the theory in orders of g and 5(= 7 ^ ~~ !)• The RG scheme will be discussed in the next 
section and details are given in the appendices. 
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III. THE RENORMALIZATION SCHEME AND THE RECURSION RELATIONS 



The renormalization group scheme we follow is based on the sine-Gordon field theory de- 



veloped by Amit et al ||16|| . The extension to the dynamics was performed by Goldschmidt 
and Schaub (GS) fl7fl . Since they presented many details of their calculations, we shall not 
repeat them here. Rather we only outline the approach and provide appendices with detailed 
explanations which complement these given by GS. 

The following renormalization constants are defined through the relations between the bare 
and the renormalized couplings: 

A) = Z D D, g = Z g g, (19) 
m 2 Q <p 2 = m 2 <p 2 R , 7o 2 2 = 7% 2 = ft = (Z$)$k> (20) 

IM> = (Z^ 1 )-^. (21) 

The subscript R labels the renormalized field variable, and the bare variable or the coupling 
constant. The terms without subscripts are renormalized constants. For convenience, we define 

h = tv- 

The renormalization of \x depends on and Zi, and no additional Z factor for renormal- 
ization is required. This is due to the FDT which implies: 

- d(t)j t (<P(x, t)(f>(0, 0)) = iA(<f>(x, t)i(0, 0)). (22) 



Here 6{t) equals to 1 as t > 0, and as t < 0. Eq. (|2l| ) is obtained by substituting Eq. (P0|) 
into Eq. (0). As we will show later, D[i will not suffer any renormalization. Therefore, it is 
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not necessary to calculate Zd in the harmonic model, which obeys the fluctuation dissipation 
theorem (FDT). 

As outlined by GS the model has a very important symmetry 4>(x,t) — > (j){x,t) + f(x), 
where f(x) is an arbitrary spatial function constant in time. As a result <p{x,t) cannot be 



renormalized and = 1 [ 17 1 to all orders in g. 



We also mention that the lattice effects and the quenched disorder in the substrate violate 



the Galilean symmetry, which provides another Ward identity z + a = 2 \TQ\ for the system 
without these effects. 

The calculations of the Z factors are exemplified in the appendices in which the explicit 
calculations of some of them are given. 
Recursion Relations : 

Once the Z-factors are known to the leading order in g, the recursion relations are obtained 
via the so-called /5-functions []18-20| : 



p> = = ^-dr )b = ( ^r }/i ' (23) 

a ( dD \ n f 9\nZ D , l 2 g^c 

Pd = <^) b = -D<-^~)b = (— A (24) 

dg ,d\nZ g 2ixg 2 

& = <-) h = -y«(-^-)6 = 25g + (25) 

,dv^ 7T7 2 9 , , 

* = ^ = -^k? a • (26) 

where subscript 6 means that all bare parameters are fixed when one performs the differentia- 
tions U-H and K is a mass scale. The renormalization of the couplings may also be related to 



the same f3 functions. Their flow under a scale change by a factor h = exp(Z) is given by minus 
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the related (3 function, in addition to the naive dependence which originates in the rescaling of 
x — > bx k — > b~ l k, and t — > b z t. The recursion relations so obtained are as follows: 

d i - °- W 

= 7T7 2 - 2 f28) 

m = |2 " 2 + £ 5^ IA (29) 

rff = - |2 - I + -5^-1"- (30) 
* = [2 -^T ]! '-(S^ 9 ' (31) 

In the next section we discusse the asymptotic scaling behaviors implied by these flow equations. 



IV. DISCUSSION : 

Recalling that the temperature of the system is T = Dp, we find that ^ = 0. Hence the 
temperature is not renormalized. The behavior of the system is governed by the renormalization 
of the coupling g. The flow of g depends crucially on the temperature. Let us define T sr = 
The recursion relation for g takes the form: 

Therefore for T > T sr g flows to zero, while for T < T sr g flows to a fixed point of order 
—5=1 — ^- with a continuous line of fixed-points [See Fig. 2]. 
We now analyze the dynamics in each phase separately. 
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A. The high temperature phases: T > T sr 

Since g — > in this phase the equilibrium properties are the same as in the high temperature 
rough phase if a surface on a smooth substrate: C(L) ~ ln(L). 

Way above T sr the mobility of the surface is finite. However, as T sr is approached the 
mobility becomes smaller and eventually vanishes at T = T sr . 

Integrating the recursion-relation we find that: 

~ H^r - 1) 2V5 (33) 

sr 

with 1J~c = 1.78. 

The dynamic exponent remains z = 2 throughout this phase, although the asymptotic 

-i 

scaling behavior is reached only on scales L > L g where L g ~ g^ 5 is the scale on which g 
decays to zero. The scale L g diverges as T — > T sr since ln(L 9 ) ~ — l) -1 . 

As long as one sits at a temperature T > T sr , the decay of g under the flow will not alter the 
asymptotic scaling behavior for L > L g , except that the amplitude in the correlation function 
depends on the bare value of g. 

In the high temperature phase, the flow of p,(l) can be calculated in terms of g(l): 

f l (l)=fi(0)e- fi So^ dl ', (34) 



where 3 = ^t^, and 



-5g(Q)e- 6 ' 
-5 — ag(0) + ag(0)e~ 



with a 



2tt 



(DA) 2 • 

Now we can obtain the macroscopic mobility, jl ~ jx{l — > oo), in the linear response regime 
(F — > 0) for the high temperature phase. 
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where § = 2^ ~ 1.78. 

As T — > T+, /i M vanishes continuously as demonstrated above. In the low-temperature 
phase g(l) flows to a finite value. Its fixed-point location changes with the temperature. More 
explicitly, the set of fixed points of different temperatures forms a fixed line in the plane of g 
and T, in which g*(T) ~ T sr — T, to first order. In this phase the scaling equilibrium properties 
and the dynamics, as well as the transport properties, are drastically modified. Most of the 
forthcoming discussions are devoted to this new, super-rough phase. 



B. The Low Temperature (T < T sr ) Super-rough Phase 

The regime of temperatures below the transition provides the most exciting new physics. 
The theoretical predictions are: (i) The correlations will change from C(L) ~ InL to C(L) ~ 
(InL) 2 . Hence the surface is even rougher than that in the rough phase at T > T sr . This 



behavior was dubbed by Toner and Divincenzo as super rough. |[21|| . They found it first in 
surface of crystals with bulk disorder, (ii) The dynamic critical exponent z now displays a 
temperature dependence. It increases continuously from its value above T sr z = 2, and z — 2 
is to first order linear in T sr — T . (iii) The linear response macroscopic mobility vanishes. The 
response becomes non linear ||, at least close to the transition, such that the average velocity v 
scales with the external force Fasti ~ F v+1 , where r] is also a temperature dependent exponent 
(related by scaling to the dynamic exponent z). 

In the following we concentrate on each of these physical manifestations separately: 

1. Super-rough Equilibrium Correlations 
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For T sr > T, g approaches a line of fixed points g* = (T sr — T)/ir (see Fig. 2). As was 



shown by Toner and DiVincenzo pi] , the correlation function of the 2-d surface acquires a 
second logarithmic factor. Since they used the static replica formalism for a model with bulk 
disorder, we explain here how this behavior is obtained within the dynamics formualtion, for 
the case of disorder in the substrate. 

The Fourier transformation of C(L,0) 2-point vertex function To t 2(q,t = 0) is defined as 
before: 

r ,2(g,* = 0) =<< h(q)h(-q) »disorder ■ (37) 

For systems with vanishing g, To^il) — ^72- Once the term u J J J dxdtdt'Vh(x,t) Vh(x,t') 
has to be accounted for, the height-height correlation function are obtained from the quadratic 
part of the Hamiltonian which contains this term. One obtains the equal time dynamic corre- 
lation, (namely the static correlation function): 

r °.* 4 = °) = + (38) 

The first term is the equal time correlation originating from the time dependent part of the 
correlation. The second term arise from the so called time-persistent part, which will be ex- 
plained in Appendix A. This can be viewed as an EA type order parameter as mentioned in 
the previous section. 

The second term is proportional to \ as the first. However, because v{l) increases with 
I (while v remains unchanged with /), it carries another scale dependence. The height-height 
correlation function, under as scale transformation by the factor b = e l , changes as: 

< Hq)h(-q) >= r 0j2 (q;v,u(0),g(0)) = e 2l T 0i2 [e l q;u,v 7 g(l)]. (39) 
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The first term in Eq. fl39|) comes from the naive dimension of < h(q)(—q) >. Using Eq. ( j38| ) 
with the renormalized values one obtains: 

< h(q)(-q) >= e 2 ^[l + ^] = -^[1 + (40) 

Apart from a finite part, v(l) increases as ~ {g*) 2 l- By choosing I such that e l q — Q — a^ 1 , 
(where Q, the momentum cutoff, can be chosen as the inverse lattice spacing), the vertex 
function, T ^(q, t = 0), is found as: 

r 0)2 (<Z,t = 0) = — 2 [A - B^\n(qa)\, (41) 
vq A v 

where A = 1, B = Consequently the static correlation function is: 



C(L) = A'\n(-) + B'\n 2 (-). (42) 
a a 



Hence, the behavior found by Toner and DiVincenzo |2T[ for the bulk disorder is also repro- 
duced in the system under consideration, where only the substrate is disordered. Experimen- 
tally, it might be difficult to distinguish hi(^) from ln 2 (^). However, the dynamical behavior in 
both phases are apparently different and their difference may be detected by the experimental 
observations. 

The behavior found above persists in the regime L « t 1 ^. 

For L » t 1 / 2 simple scaling implies a dependence of (lnr) 2 on r. The derivation of the 
intermediate behavior of C(L,t) is beyond the scope of this paper (it requires the knowledge 
of To^{q, w). However, simple physical considerations hint very strongly that its behavior is of 
the form C(L,t) ~ {\nLf[r/L z ]} 2 . 
2. The Dynamic Exponent 
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To calculate the value of the dynamic exponent z in the low temperature phase one should 
look at the recursion relations of D and Jx. To locate the fixed point one can require that both of 
^ and ^ to be equal to zero. Then z is obtained as z = 2 + Ay/c\8\, where the fixed point value 
of g, g* — was inserted. As far as we know, this is the first example in which the present 
dynamic exponent z varies with temperature continuously (besides the random anisotropy XY 
model, which is described by the same theory studied by GB). The physical implications of the 
increasing z can be understood from the fact that the relaxation time to reach the equilibrium 
state is longer as the temperature is lowered below T sr . This is to be expected in the phase 
where the disorder is relevant, and the surface turns super- rough as explained in the previous 
section. Since the surface stretches itself to find the configurations with the lower free energy, 
it will leave these locations slowly. The slower dynamics implies also a graduated increase in 
the averaged free energy F(L) barriers, associated with a scale L, due to the disorder. 
3. The Non-Linear Response 

In the last section, addressing the high temperature rough phase, we have found that the 
linear response mobility vanishes as (T — Tsr) 1 - 78 , when T is reduced to T sr . Below T sr the 
linear response mobility p,M vanishes. 

As we show in the following the response becomes non-linear. Again the physical origin of 
this behavior takes roots in the preferred configurations of the surface which are local minima 
of its free energy. Applying a small force F will not move the surface in a uniform velocity. Ac- 
tually a somewhat similar situation occurs in the smooth phase of a surface growing upon a flat 
substrate, where the mobility jump from a finite value to zero at the roughening temperature. 

If the pinning is even stronger the surface grows by activation of higher islands. If the 
substrate is disordered, the preferred and inhomogeneous locations of the surface are enough 
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to slow the motion and then to cause the linear response mobility to vanish. However, it still 
allows for a uniform motion with average velocity which vanishes as F 1+r > (rj > 0) when F — > 0. 
Hence, the velocity vanishes faster than F. The force F is a relevant field which increases as 
F(L) ~ FqL 2 with the length scale. Even for small Fq there is a scale L* ~ aF for which 
the scaled force is of order 1. Namely it is not a negligible quantity. On the scale, L > L* the 
behavior is not "critical". The large force moves the surface with a uniform velocity. The ratio 
between the force and the velocity is determined by the mobility \i at the boundary L* between 
the "critical" (L < L*) and the "non-critical "(L > L*) regimes. Within the scaling picture L* 
can serve as a "cutoff length ". Hence it is the mobility l^(L*) for a piece of the interface with 
linear extent L* that is determining the mobility of the whole surface on scales L > L*. Note 
that fi(L) is not exactly p,(L*) but the ratio between them is finite since no " critical scaling " 
holds for L > L*. The scaling of p,(L) for L < L* may be derived from its definition as: 

Under rescaling L — > L/b we already know that t — > t/b z and F — > Fb 2 while /i does not 
suffer any renormalization. To make both sides scale similarly we must have /i — > fib 2 ~ z and 
defining /i ~ /U &~ 2t ' we obtain: 

77=(^-2)/2 = -2 > /o5 = 1.78|(J| (44) 

That implies that fi(L) on a scale L < L* is made smaller by (^)~ 2ri with respect to its bare 
value. In particular we find: 

= Vo(~r 2v = VoFS = VoFo /2 ~\ (45) 
a 

where we have used relation between L* and F . 
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So we have identified the dependence of \lm on F from which we obtain the averaged 
velocity: 

v ~ fi F^F = fi F^ = fx F* /2 . (46) 

Thus as temperature is lowered the velocity (for the same tiny force Fq) becomes smaller. 
How far below T sr these relations hold? The scaling picture is based on a local equation of 
motion of h. Therefore, it implicitly assumes the existence of single solution for the equation 
of motion in the limit of vanishing uniform force. Just below T sr this is a valid assumption 
since even if more than one minima exists the scale associated with the difference between the 
minimizing configurations (which diverges as T —>■ T sr ) is larger than the scale we discussed 
here. At lower temperature this may no longer hold. The existence of multiple minima might 
be felt on the relevant scale (i.e. L* ). In the regime where many minima are relevant the 
dynamics will be activated. In other words, the slower processes are going to be related to 



the height of the barriers between these minima. Some works p2| , pT3| have been devoted to 
estimating these barriers and drawing the conclusions based on " activated dynamics". It is 
not clear, however, how reliable these heuristic estimates are (e.g. the barriers are identified with 
the fluctuation in the minima of the free energy.) Unbounded barriers between configuration 
unrelated by symmetry will also lead to broken ergodity which can be reflected in "replica 



symmetry breaking" . Such a possibility was found recently within a variational approach [^3 . 



C. T at and just below T sr 

The discussion we presented so far for T < T sr applies on length scales for which g is already 
close to its fixed-point value g*. If the temperature is very close to T sr this scales became very 
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large and it will be necessary to account for the crossover regime. The recursion relations can 
be integrated as before and yield for the scale dependent mobility: 

HI) = mil + ^(0)^[(^) 2151 - l]}- 2 ^. (47) 

^From this expression we see that if (^) 2 ' 5 ' >> 1 the results of the previous section apply. If 
the force F is not small enough for the associated scale L* to satisfy this relation the condition 
between the velocity and the applied force changes to: 

^-{l + ^^^'-l)}- 1 - 78 (48) 
At T = T sr (5 = 0) this expression yields: 

i = {1 + *m HF)) ->.» (49) 

Namely v/F vanishes as F — > due to logarithmic corrections. They originate from the effect 
of g which is marginally irrelevant and decays to zero so slowly that it still causes v/F to vanish. 

V. SUMMARY 

To summarize the main conclusions of our investigation on the effect of disorder in the 
substrate on the surface dynamics: There exists a phase transition between a high temperature, 
rough, phase, and a low temperature, super-rough, phase. 

In the rough phase correlations and response functions have the same scaling properties 
as in the pure case. The disorder and the periodic potential are irrelevant in this phase. As 
the transition temperature T sr is approached from above the macroscopic mobility vanishes 
continuously as (T — T sr ) 178 (contrary to the flat substrate in which it has a finite jump). 
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The properties of the low temperature phase are unique. The height-height correlations are 
C{L,t) ~ [InL] 2 as L « r 1 ^, and C(L,t) ~ [lnr] 2 as L » r x l z with z = 2 + 4^(1 - £-) 
to first order in (1 — T/T sr ). Namely the dynamic exponent increase continuously from its 
Gaussian value of 2 as the temperature is lowered below the transition. 

The linear mobility has a scale dependence which causes it to vanish on large scale. We 
have shown that this scale dependence results in a non-linear relation between the applied force 
and the average velocity: v ~ F 1+v , where r\ = 1.78(1 — is also a temperature-dependent 
exponent. 

All these results also apply to the 2D vortex-glass system in a film of type-II superconductors. 



We discussed these implication elsewhere []14 |. The most important one is the non-linear 



relation between the voltage V and the current I: V ~ I 1+v with 77 given above. 

Our RG calculations give the same static behavior obtained by the replica approach pre- 
senting the unbroken symmetry between replica. 

Other works have shown the symmetry to be broken within a non-perturbative variational 
gaussian approach which is equivalent to the N — > 00 limit. It is an unsettled matter whether 
the symmetry is indeed broken for N = 1. 

Preliminary numerical results for a 2D vortex-glass [27] and a random growth model 



show a transition in the dynamic properties at a temperature within 10 % of the analytic RG 
result. Looking at the static correlations below T sr , however, the possibility of the replica 
symmetry breaking may not be excluded. However, It is too early to draw firm conclusions 
from those preliminary results. It is to be expected, therefore, that more analytical and nu- 
merical works would be necessary to reconcile the different results and to reach a complete 
understanding of this exciting and challenging problem. 
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APPENDIX A: RESPONSE AND CORRELATION FUNCTIONS 

Two kinds of propagators, the response function and the correlation function, arise from the 
free (Gaussian) portion in the effective action in Eq. (|TTD (consisting of a quadratic form in the 
field (f) and the auxiliary field <fi). One can calculate the free response and correlation functions 
directly from the free part of the action Sq in Eq. (|I~5"D. In the momentum and frequency 
representation, they are: 

(<f>(q, uO0(-£ -u)) = ,2*^. , (Al) 

jj,{q z + m z ) + iuj 

(<f>{q, -U))) = 2 , 2D 2M2 , 2 ' ( A2 ) 

where m is the mass of the field <fi. We have introduced the mass m to regularize the in- 
frared divergences which appear in the upcoming loop calculations. In the momentum and time 
representation, they are given by: 

(mtm-q,t')) = 9(t-t')e-^ 2+m ^ t - t '\ (A3) 
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<<Kg,t)<K-g,t')> = -2^-e-^ 2 )\^\ (A4) 
qz _|_ m z 

In this manner, both the free response and free correlation functions possess a mass depen- 
dence in their denominators as one can see from Eqs. (|^T|)-Eqs. ( |A4[ ). 

The short distance cut-off is introduced for regularizing the ultraviolet divergence in 2 
dimensions. The regularized Co and Rq are: 

c ^ ') = / Wr^'^A^'^ < A5 > 

^^'^ = " (A6) 

Furthermore, the following equations are very valuable to extract the asymptotic behavior 



of free propagators \ \L% : 



1 1 1 x 2 + a 2 

C (x,t) ~ ln(mlX Q \t\) C — — , as m 2 (x 2 + a 2 ) « m 2 r.\\t\ « 1 (A7) 

An An An AX\t\ 

C (x = 0,t) = -^\n(2V~cm 2 \\t\) + 0(-^-), (A8) 

An Ao 1 1 1 

where C (x,t) = DfiC (x,t), and C is the Euler's constant ~ 0.5772. In the limit moA|i| << 
m 2 (x 2 + a 2 ) << 1, the equal time correlation behaves as: 

C (x, t = 0) = -— \n(cm 2 Jx 2 + a 2 )) + 0(x 2 ), (A9) 
An 

The zero order fluctuation-dissipation theorem (FDT) relates the response function to the 
correlation function as follows: 

1 9(t)^C (x,t) = R (x,t). (A10) 



fi 2 D y ' dt 
1. Basic diagrams 
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In Fig. 3 the wavy line represents the auxiliary field 0; the straight line represents the field 
(f). The dash line is set to separate two different time coordinates. The dot points represent the 
abbreviation of the other lines. 

2. Free propagators: 

Correlation function Co(x,t) and response function Ro(x,t) are shown in Fig. 4. 

In general, the equal time correlation function, C(x,t = 0) =< <f)(x,t)<f>(0,t) >, is identical 
to the static correlation function < <f>(x)<f>(0) > averaged by the Boltzmann weight. As a special 
case for the connection, one can refer to Eq. ( |A4j) . On the other hand, the response function 
is related to the reaction of the system to the external probe, say P(x,t). The perturbed 
Hamiltonian will result in TC — J d d xP(x,t)<f)(x,t), and the additional term in the Hamiltonian 
is tantamount to adding P(x, t)no4>(x, t) in the MSR action. In other words, the conjugate probe 
P(x,t)no will couple to the response field (f>(x,t). Intuitively, the linear response function can 
be defined as the ratio of the strength of the reactive field to that of the probe: 

Rp hy {x-x ,t-t) = ^ =^ Jjn, J\ = M<P{x,t)(f){x ,t)). (All) 

where fj,oP(x, t) = J(x, t) (note that R p h y is different from R by a factor /i and the latter 
is used in the calculation here). This arises from the various representations of 0, which is 
associated with the unphysical degrees of freedom in multiplying the Langevin equation by an 
arbitrary constant. 

Consistently, the response function will manifest causality. There is no response before the 
external probe is applied, i.e., R(x,t) = for t < 0. 

The relation between the static and the dynamic properties are through the correlation 
functions. In the random systems one distinguishes between equal-time correlations which 
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correspond to the static correlations: 

< (j)(x,t)(/)(x ,t) >= [< (j){x)(j){x) >}av (A12) 

where [..]av stands for the average over disorder. In the replica calculation it is given by the 
diagonal term < (j) a (x)(f> a (x) >. Time persistent part of the correlation which correspond 



to the Edwards- Anderson (EA) type of static correlation [|25|j : 



[< 00, t)(f)(x',t ^'^Jav = Qea{x -x), (A13) 

In the replica language, if there is no replica-symmetry breaking, it is given by the non- 
diagonal part: lim„_ >0 [< 4> a {x)(j) l3 {x) > a ^f3\AV- 

APPENDIX B: PERTURB ATIVE EXPANSIONS USING FEYNMAN DIAGRAMS 

To facilitate the RG calculation, one can introduce lPI(l-particle- irreducible) vertex gen- 
erating functional T[M]. In analogy to the static case, one can define the vertex functions 
through T. For example, the vertex function r^i is define by gjj^j ||20|| . Here the calculation 
is focused on expanding the bare parameters. On the other hand, one also can expand the 
theory in terms of renormalized parameters. Up to 2nd-order perturbation, the 1PI vertex 



generating functional are expanded as [26 



T[M] = ^(M,G M)+P[M], (Bl) 



P[M] = (V[0 + M]) - \{{V 2 [<P + M]) - (V[0 + M] 2 >) } - {( ^+ M] )o, 

^0(^4^)0} + 0(V% (B2) 
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where G$ 1 is symbolizing the 2x2 free propagators matrix, composed by Rq and Cq, and the 
bracket, <>o, stands for the average taken with respect to the free action. The second term 
in the right hand side of the equation is the sum of the connected diagrams up to the order of 
V 2 , and the third term is serving to remove all connected diagrams in the second term which 
are not lPIs. 

To proceed with the calculations systemically, it is convenient to introduce the vertex func- 
tions here. The bare and renormalized vertex functions can be related by factors of Z^, Zg. 
For instance, 



r iv,i,(9> w; X R , m R , k) = (Z^) 2 (Z$) 2 T NjL (q, u; A , m , a), (B3) 

where A# and Ao label renormalized parameters (gn,[i, ■ ■ •) and bare parameters (go, no, ■ ■ •), 
respectively, q and u are the external momentum and frequency, respectively. In the corre- 
sponding vertex function, a is a short-distance cutoff, and k is a mass scale. Here Tn,l stands 
for the vertex function with L external <ft lines and N external </> lines. The factors, Z^ and Zi, 
are set to remove the divergent parts of the vertex function V. 

APPENDIX C: THE CALCULATION OF Z- A 

<p 

Prior to the calculation, we make some remarks about the notations of the parameters. 
Since mR and 7^ satisfy Eq. ( |19|) and Eq. (^) with Z^ being 1, we have m = m and 70 = 
7. Thus there should be no confusion if we use them interchangeably. Before entering the 
calculations of the Z factors, one should make reference to the relations between the various 
vertex functions. For the calculation of the Zi factor, one should focus on the vertex function 



17| ) Ti i. Inferred from Eq. ( |B3| ) , their renormalized and bared counterparts are related by: 
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rji(ga;; fi,v, g) = Z^T hl (quj; n ,m ,go) = Z<f,T ltl (quj; fiZ- 1 , m, Z g g). (CI) 

The contribution from the action of the first order in g to Z^ can be calculated by considering: 

5 2 

5M(x, t)8M(6, 0) 

where 



(V g [<f> + M,<f> + M]), (C2) 



2 2 

Vg[ij>, <p] = J J d 2 xdtdt'0(x, t)$(Z, t') cos[7(0(x, t) - (j)(x, if))]. (C3) 

Fig. 5. shows one representative of the diagrams corresponding to T^i in Eq. (|U4j) . The 
dots inside the circle stand for the contractions of <f>(x, t)s and <f>(x, 0)s, which results in e' y2Co ^ 0,t ^ 
in Eq. (|C4]). The dots near the left "ears" represent the contractions of 4>(x,t)s and (p(x,t)s, 
which result in e~ /y2c °^ ' ' > in Eq. ( |C4| ) . The roles of the dots near the right "ears" are similar. 

In the frequency representation, the integral coming from the contractions of inner lines can 
be written as: 



/oo 
e iuJt dtR(0,t)e 
-oo 



-7 2 C (0,0)+7 2 Co(0,i) 

' — oo 



Jo 



c/te^e" 7 Co(0 ' 0) } (C4) 



e - 7 2 co(o, 0) r dte ^t d ±d 2 2Com _ 

Jo fi z D dt 

e -^c o{ o,o) luJ f 00 dte i^c (o,t) _ i]( J_) (C5) 
Jo fi 2 D 

2„2\6 



Combing with the prefactors and using g = go(cm ra f) , we have 



■ \ 2 2 H'O UL 

cm 



fx D Ja 2 /x 2V cm 2 fj, t 



(C6) 



Finally, Z^ is found as: 
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^= 1 + 2TS ln < KV >> < C7 > 

where k = cm 2 is the scale, at which one can impose the prescription for the vertex functions 



corresponding to the renormalized parameters [26 



APPENDIX D: THE CALCULATION OF Z D 

To calculate the factor Z^, we consider the vertex function r 2j o- The bare vertex function 
r 2 ,o is related to r^ through Zi. They can easily be seen from Eq. (|B3|) : 

rj (£ w; g) = zlr%{q, u- z^, z g g). (di) 

The associated diagram is illustrated in Fig. 6. After taking the contraction of the inner 
lines, we are left with: 

~ 9J ^j f M(z,t')M{x,t), (D2) 

times 2 < cos{7[0(x,t) — <f>(x, t )]} >o, where the contents inside the bracket are averaged with 
respect to the free action. The latter term contributes to the renormalization of two-point 
vertex depicted in Fig. 6. The calculation of this term can be performed with ease as shown 
below. In the time representation, it reads: 

_ 2 9^ e -^c (o,o) f°° dt(f ?Com e M (D3) 
2a 2 J-oc 

In the limit — ^ 0, it reduces to: 

_ 2 £^7o e - 7 »cb(o,o) r dte ^c (o,t)_ (D4) 

2(2 J — oo 

With the aid of the formula given in Appendix A, the divergent part of this term can be 
obtained by proceeding in the same manner as in the previous section: 
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_ 2 g^7| e - 7 2 Co( o,o) 2 ln(cm 2 Q2) + finite termg _ (D5) 
2a 2 2^cm 2 [i 

Then, 

-2/V = (^) 2 {-2 J D ^ + 7 2 ^/ioln(cm 2 a 2 )}. (D6) 

Finally we inherit Z D = 1 + 7 2 ^f ln(cm 2 a 2 ). There should be no confusion if we put either the 

bare or the renormalized values of D and \i in the denominator in the above equation. Up to 

this order, there will be no difference between these two possibilities. 

APPENDIX E: THE EXPRESSION OF Z g AND THE RENORMALIZATION OF v 

It is a tedious calculation to find Z g by considering the vertex function T 2! o{q,uj). Since 

quenched disorder is present in the system, one may instead consider r 2j o(<f, t\\ —q,t 2 ), in the 

limit \ti — t 2 \ — > oo. 

To obtain the scaling equation of u, one should consider the q dependent part of 
^2,o(q, ti\ —q,t 2 ) with \ti — t 2 \ — > oo. The contractions of fields which connect the points t\ and 
t 2 , make no contributions. The q dependent part of r 2i o up to order g 2 involves 7 diagrams, 



which are classified into 3 sets as illustrated in GS | 17| . Those terms can be expressed as the 



first seven terms in (7.3) of GS [17]. The first set (referred to figl. (a) in GS) contains only one 
diagram as shown in Fig 7 (A). The second set contains three diagrams, corresponding to part 
(B), (C), and (D) in Fig. 7. The third set can be obtained from the second set by swapping 
(xi,ti) < > (x 2 ,t 2 ) and (xi,t") < > (x 2 ,t). 



In the following, we only present the calculations which are not explictly given in GS [fPT|. 
In the calculation, we only concentrate on the time dependent part, neglecting the prefactors. 
This should not be a cause for confusion when one retrieves them later. By changing the 
variables t' = t 1 — t' and r" = t 2 — t" the first term is recast into: 
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/oo f roc a i ii i ii 

dr / dr R (x,t )R (x,t ) sinli7 2 [— C (x, r ) — C (x, r )], (D7) 
-oo J —oo 

where x = x\ — x-i- 

With the identity of FDT for R$ and Co and the integration by parts, it is simplfied into: 
-1 r 00 r 00 , , , rr d . , „ „ , d 



(jj?D~i 2 ) 2 Jo Jo 



f f rfr'rfr''{[— sinh7 2 (7o(a;,T')][— 7Cosh7 2 C7o(a;,r'')] 
Jo Jo dr dr 



+ [ cosh7 2 C* (a;, r')] [^77 sinh7 2 C (£, r")]} 

~~ sinh(7 2 Co(x,0))[cosh(7 2 C (x,0)) - 1]. (D8) 



By changing the variables r = t\— t and r" = t' — t" in the second term (see part (B) in 
Fig. 7 ) and r = t" — £ and r" = t 2 — t" in 5th term, they become: 

/OO /"CO 
dr / dr R (x,r )R (x,r ) exp[7 2 C (0, r +r ')]sinh7 2 [C (x,r')-^o(x,r")]. (D9) 
-00 J —00 

By changing the variables r' = t\ — t' and r" — t — t" in the third term (see part (C) in 
Fig. 7) and r' = ti — t" and r" = t" — t' in the 6th term, they become: 

/CO r OO 

dr' / dr"R (x, t')R (0, t + r") exp[7 2 C (0, r + r ')] cosh7 2 [C (x, r') - <7 (a;, r ')]. 
-00 J — oo 

(D10) 

By changing the variables r' = t ' — t' and r" = t — 1\ in the 4th term (see part (D) in Fig. 
7) and r —t' — t" and r" = t — in the 7th term, they become: 

/OO . /*0O 11 1 1 11 ill l 11 

dr / dr i?o(^, t )Rq(0, — r — r ) exp[7 2 C (0, r +r )] cosh7 2 [C (x, r ) — C (x, r )]. 
-00 J —00 

(Dll) 



With the FDT identity, Eq. ( |D9| ) can be futher simplfied into: 

W^Y So dT So dT " ex P^ 2 ^o(0,r' + r'')][^cosh7 2 C'o(a;,r')^sinh7 2 C , o(a;,r'') 



(/i 2 J D 7 2 



77 cosli7 2 Co(x, r )— - sinh7 2 Co(x, r )]. (D12) 



dr" "' dr 
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By using the FDT and the integration by parts, one can reduce Eq. ( P10| ) into: 

J dr ' J , dr " exp[7 2 C (0, r' + r'')][^ 7 cosh7 2 C' (x,r')^ 7 sinh7 2 C'o(a;,r'') 



- -^77 cosh7 2 C (a;, r")-— sinh7 2 C (x, r)\ 
dr dr 

+ the boundary term . (D13) 



Simliarly, Eq. ( [L> 1 1| ) is recast into: 

1 r°° , / r~ T , // r „ ,„ / »„ r d , 9 _ , d 



(/i 2 Z>Y 2 ) 2 Jo 



/ dr' / dr" exp[7 2 C (0, r' + r")][— cosh7 2 C (x, r')— 7 sinh7 2 C (a;, r") 
Jo ar ar 



— cosh7 2 Co(x, r )-— sinh7 2 Co(x, r )] 



dr ar 
+ the boundary term . (D14) 



Combing Eq. (|D12|) , Eq. ( P13| ) and Eq. ( P14|) terms excluding the boudary terms, we have: 



T^Y Jo dT ' I dT " ' ex P^ 2 ^o(0,r' + r'')][^cosh7 2 C (^,r')-^sinh7 2 C , (a;,r''' 



- -^ 77 cosh7 2 C (a;,r")-^ 7 sinh7 2 C'o(a;,r')]. (D15) 
dr dr 



The above integral vanishes since the integrand is antisymmetric with respect to the trans- 
formations t — > — r" and r" — > —r. 

The boundary terms in Eq. ( |D10| ) and in Eq. ( [01 1| ) cancell each other also. For example, 



the boundary term in Eq. (|D10| ) reads as: 
, d 



roc rl roo 

/ dr— [sinh 7 2 Co(x,r)]+ / dr exp[ 7 2 C (0, 0)] 
J dr Jo 



[— - sinh7 2 C (x, r ) cosh7 2 C (x, — r ) — - cosli7 2 (x, r ) sinh7 2 C (x, — r )]. (D16) 

ar ar 



To sum up, the only contribution is from Eq. ( |D8| ) . Extracting the singular contribution 



I7yi9|] gives v h = v + s^tSf ln ( K2a2 )- 
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Figures Caption 



Fig- 1 

A two-dimensional cut (along a lattice plane perpendicular to the disordered 
substrate) of the three-dimensional system. 
Fig. 2 

The flow diagram for g(l) with two different fixed-lines for T larger and smaller than T s . 
Fig. 3 

The Feynman diagram representing <p(x, t)(f>(x, t') cos(0(x, t) — <p{x, t')). 
Fig. 4 

The lines representing the free correlation and response functions. 
Fig. 5 

The Feynman diagram for the vertex function 
Fig- 6 

The Feynman diagram for the vertex function J^o- 
Fig- 7 

4 Feynman diagrams contributing to the renormalization of v. 
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